Fluctuation-induced pair density wave in itinerant ferromagnets 
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Magnetic fluctuations near to quantum criticality can have profound effects. They lead to char- 
acteristic scaling at high temperature which may ultimately give way to a reconstruction of the 
phase diagram and the formation of new phases at low temperatures. The ferromagnet UGe2 is 
unstable to p-wave superconducting order - an effect presaged by the superfluidity in He 3 - whereas 
in CeFePO fluctuations drive the formation of spiral magnetic order. Here we develop a general 
quantum order-by-disorder description of these systems that encompasses both of these instabilities 
within a unified framework. This allows us to demonstrate that in fact these instabilities intertwine 
to form a new phase, a pair density wave. 

PACS numbers: 74.40.Kb,75.30.Kz,71.10.-w,74.20.-z 



INTRODUCTION 

The physics of metallic systems tuned near to contin- 
uous, zero-temperature phase transitions has provided 
a rich vein of experimental surprises and theoretical in- 
sights [IH3] ■ Fluctuations in such systems drive new crit- 
ical behavior that inextricably involves quantum dynam- 
ics, described by Hertz [4] as quantum criticality. Over 
the past decade, it has become apparent that quantum 
fluctuations have an even more profound effect on the 
magnetic phase diagram [5HTU]. A breakdown of the 
standard Moriya-Hertz-Millis gj El IE] model of itin- 
erant magnetic criticality occurs as the quantum crit- 
ical point becomes occluded by new fluctuation-driven 
phases. Magnetic fluctuations transverse to the polar- 
ization drive many materials to display a first order 
ferromagnetic transition [13] . spatially modulated mag- 
netic order - for example CeFePO [14], NbFe2 [15] . 
Yb(Rh . 73 Co . 27 ) 2 Si 2 [16 and possibly Sr 3 Ru 2 7 [TT] 
118] , or higher angular momentum particle-hole pair- 
ing PjHrjO]- Moreover, in materials such as UGe 2 pTrE5] . 
URhGe [H], and UCoGe [25J E§] longitudinal magnetic 
fluctuations can drive p-wave Cooper pairing in a similar 
manner to He 3 superfluidity [27H29] . 

In this Letter we study the minimal model of freely 
dispersing electrons interacting through a contact repul- 
sion. We develop a general quantum order-by-disorder 
approach to calculate the effects of fluctuations on the 
free energy, leading to the phase diagram in Fig.[l] Near 
to the itinerant ferromagnetic quantum critical point the 
electron gas forms a previously unknown phase, a pair 
density wave where the magnetic order forms a spiral 
modulation coexistent with a Larkin-Ovchinnikov-Fuldc- 
Ferrell-like spatial modulation of the d- vector of a p-wave 
superconductor [30H32] , 

The essence of our order-by-disorder approach is to 
expand the free energy in fluctuations about some pu- 
tative order. The order modifies the electron disper- 
sion and so the energy of the quantum fluctuations. 
This offers an alternative, physical view of the non- 
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FIG. 1. (Color online) The phase diagram of Hamiltonian, 
Eqn.[I]as a function of temperature and interaction strength. 
Our new phase, the pair density wave, is highlighted by the 
red text, g is the interaction strength, u-p is the density of 
states at the Fermi surface, T temperature, and Tp the Fermi 
temperature. The magnetic background to this phase diagram 
is as determined in Ref. [35], allowing for a resummation of 
the leading divergent contributions to the free energy. At high 
temperatures, there is a continuous transition from ferromag- 
net to paramagnet. The temperature of this reduces with 
decreasing interaction strength until the red tricritical point. 
At this point, there is a Lifshitz transition from the ferromag- 
net into a magnetic spiral (green dotted line). Finally, the 
spiral magnetic order gives way to paramagnetism in a first 
order transition. Superconductivity in the absence of magnetic 
order occurs below the blue dotted line. This transition tem- 
perature rises dramatically in the presence of magnetic order 
due to mode-mode coupling. Where superconductivity over- 
laps the spiral magnetic order a Larkin-Ovchinnikov-Fuldc- 
Ferrell-like pair density wave state is formed, highlighted by 
the red text. 

analytic corrections that overturn Moriya-Hertz-Millis 
theory leading to the opportunity to undergo a first or- 
der transition [5J [5j H3J 1231 131] > spatially modulated mag- 
netism [35] [35], a nematic phase [5U], and p-wave su- 
perconductivity [37MT] . Near to the quantum critical 
point the fluctuations self-consistently stabilize the or- 
dered state. This mechanism is similar in spirit, but sub- 
tly different in details, to several others: the Coleman- 
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Weinberg mechanism of field theory [32] , quantum order- 
by-disorder [43], and the Casimir effect [44]. From the 
latter perspective, the Fermi surface may be viewed as 
providing boundary conditions for quantum fluctuations 
through Pauli exclusion. Modifying the Fermi surface 
by the introduction of order modifies the spectrum of 
the fluctuations. This perspective emphasizes the central 
role played by the underlying fermionic statistics (which 
is absent in the conventional application of these ideas). 
However, we have chosen to use the term quantum order- 
by-disorder to be closest in spirit to the conventional us- 
age. 

In this Letter we first revisit the quantum-by-disorder 
approach demonstrating how transverse magnetic fluctu- 
ations can generate a spatially modulated magnetic phase 
and extend the analytical formalism to apply near zero 
temperature. Secondly, we introduce a superconducting 
instability driven by longitudinal spin fluctuations into 
the formalism and elucidate the phase diagram. Finally, 
we discuss the consequences of the intertwined spiral and 
superconducting instabilities for experimental systems. 



QUANTUM ORDER-BY-DISORDER IN THE 
ITINERANT FERROMAGNET 



We first review the order-by-disorder approach to ad- 
dress the magnetic phase diagram ignoring any possi- 
ble superconducting instabilities. Our starting point is a 
Hamiltonian of freely dispersing electrons with = k 2 /2 
interacting through a contact repulsion g reflecting the 
screened Coulomb interaction present in metals 



ergy is given by T = J" MF + -^fluct with 
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(1) 



Here n£ is the fermion occupancy in momentum space, 
n° in real space, spin a G {t,l}, an d we adopt atomic 
units with H = m = 1 throughout. This model shows 
a continuous ferromagnetic transition in mean-field, the 
zero-temperature terminus of which Hertz termed the 
quantum critical point [4]. Hertz realized that quan- 
tum fluctuations could modify the critical behavior at 
this quantum critical point, however it has recently been 
shown that fluctuations can have an even more dra- 
matic effect, reconstructing the phase diagram [5] 0- 

Eol [331 ESI EH] . 

The analysis of these effects proceeds as follows: first 
we consider a background spiral magnetization, M q = 
M(cos q • r, sin q • r, 0) that modifies the mean-field elec- 
tron dispersion e£(M q ) = e k + er\/(k • q) 2 + (gM) 2 . The 
mean-field and fluctuation corrections to the free energy 
are then calculated in the presence of this background 
spiral. To leading order in the fluctuations, the free en- 



Jmf = -T]Tln[c^( M o>-^/ T + 1] + 9 -M 2 , (2) 
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where /i is the chemical potential, T temperature, and 
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The expectation is taken in the thermal ensemble and 
the states |s) are virtual intermediate states generated 
by the action of the Hamiltonian. The fluctuation con- 
tribution to the free energy is then calculated following 
the prescription of Refs. [33j [33] , with the result [45] 
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Ref. |35| evaluated the J-auct numerically at each value 
of M q . This revealed that the low-energy phase space 
of fluctuations is enhanced by uniform or spiral mag- 
netic order. Good agreement with variational Monte 
Carlo calculations [35] encourages us to develop a more 
amenable Ginzburg-Landau expansion near to the tri- 
critical point [3D]. The free energy is a functional of 
the mean-field dispersion in the presence of the modu- 
lated order parameter, J"(M q ) = J"[e k (M q )] = (T[M 2 + 
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(T[9 2 q 2 ]) where 6 = k • q/|k||q| and 



) indi- 
cates an angular average. In a Taylor expansion of the 
free energy in M 2 and q 2 , terms are then linked by geo- 
metric factors stemming from the angular average 

T= ^+^/3g 2 + ^ 7 ^Af 2 + (/3+ 7(? 2 )A/ 4 + 7 M 6 , (5) 

so the expansion contains only three coefficients a, /?, and 
7 that are functions of the interaction strength g and the 
temperature T . Reconstruction of the phase diagram is 
driven by a logT contribution to the M 4 coefficient j3. 
The analytical form of these parameters is known near 
to the tricritical point |20j . 

We now seek a superconducting instability that occurs 
at very low temperatures away from the tricritical point. 
In order to access this, we must extend the region of va- 
lidity of our expansion to lower temperatures away from 
the vicinity of the tricritical point. The leading fluc- 
tuation contribution, M 4 In T is accompanied by higher 
order terms proportional to M 6 /T 2 , M 8 /T 4 etc. These 
occur with pre-factors allowing them to be re-summed 
exactly to give the following contribution to the free en- 
ergy 
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where v-p is the density of states at the Fermi surface 
and C is the Euler constant. This is consistent with 
that found in previous studies [35]. In practice, 
sub-leading terms that are less divergent in T have an 
important effect upon the position of the phase bound- 
aries at low temperatures away from the tricritical point. 
We find that they enhance the region of the phase dia- 
gram occupied by the spiral phase making it consistent 
with the broader region seen in variation Monte Carlo 
calculations [35] - Unfortunately, re-summation of these 
sub-leading divergences does not yield a simple analytical 
form and we have not included them here. 

SUPERCONDUCTIVITY IN QUANTUM 
ORDER-BY-DISORDER 

We now wish to incorporate the possibility of Cooper 
pairing. We focus upon Al type pairing of spin-up elec- 
trons into a p-wave state [46], as it is the most favor- 
able given the background magnetic order. Other pairing 
symmetries can be treated similarly. 

The superconducting order overlaps with the 
fluctuation-induced spiral modulation so we must 
construct superconducting order that can be treated 
self-consistently in the order-by-disorder scheme. The 
natural choice is a p-wave state in the basis that is 
diagonal in the presence of the background spiral, i.e. 
in the basis where electron spins are polarized locally 
parallel or anti-parallel with the background spiral 
order. This leads to a state in which the d-vector 
of the p-wave superconductor rotates in concert with 
the magnetization. The state may be referred to as 
a pair density wave as it involves modulation of both 
magnetic and superconducting order, or alternatively 
Larkin-Ovchinnikov-Fulde-Ferrell-like since it involves 
a modulation of the superconducting order parameter, 



but in its SU(2) rather than U(l) sector. 

As it stands, the quantum order by disorder procedure 
requires that the emergent order changes the mean-field 
dispersion. In the case of a contact repulsion, supercon- 
ductivity does not have this effect. One way around this 
(used to study the possibility of nematic order in [5D]) is 
to include a field conjugate to the superconducting order. 
After calculating the appropriate free energy and its fluc- 
tuation corrections, a Legendre transformation recovers 
the Gibbs free energy. Here, we use an equivalent self- 
consistent approach. The general scheme is as follows: i. 
after having first diagonalized the electron states in the 
magnetic background, we add and subtract a term 

E( Ak5 -kt 2 k t + c - c ) (7) 

k 

in the free energy, where cj^ creates a state with spin po- 
larized parallel to the background magnetic order, which 
may be spiral in certain regions of the phase diagram, ii. 
The terms Ho — J2k (A^c 1 ^ + c.c) can be diagonalized 
by a Bogoliubov transformation, iii. The contributions 
of the terms "Hint + 5Z k (A k ctct + c.c) can be calculated 
perturbatively in the Bogoliubov basis. For BCS the- 
ory, the resulting free energy has the same gap equation 
as the conventional treatment. As we now show, using 
this approach in the case of a contact repulsion recovers 
spin- fluctuation mediated pairing |37H41] . 

The superconducting order does not change the expec- 
tation of the contact interaction and so only enters the 
mean-field free energy through the kinetic term. Fluctu- 
ation corrections are obtained following the prescription 
developed for the fluctuation driven spin spiral; the fluc- 
tuation contributions to the free energy are calculated 
from Eq. ^ , though now in the diagonal Bogoliubov ba- 
sis. Overall, after expanding to order |A| 2 , we obtain the 
free energy [47] 



J-(M,A) = -T^ln(e-WT + 1) + ^!_ 29 2 £ 

k,CT k 1+ k 2 £ ki + e k 2 e k 3 £ k 4 

=k 3 +k 4 

-£^(24-l)V£|^(l-2n k ^)|^^ (8) 

k Z ?k,t k q ^k+q,t ^?k,f k z ?k,t 

where £ k(J = e k — ji. x an( l S are the susceptibility and self-energy evaluated at finite M q and spiral vector q, and 
are given by 

p £ p+q £ P U p : q £ k + £ p-q £ k-q £ P 



I 

The first line of Eqn. [8] that is independent of A is pre- cisely the free energy of the ferromagnet derived earlier. 
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At finite A, the first term of the second line contains the 
mean field response of the Fermi surface to Cooper pair- 
ing, and the subsequent terms arise from the fluctuation 
corrections. 

The gap equation resulting from this is given by ex- 
actly the same expression as the Eliashberg equations 
obtained in spin-fluctuations theory. The resulting tran- 
sition temperature for the pair density wave can be writ- 
ten as |BTH4T] 



T c = exp 

7T 



(l-0 e reeSt( k , ek ))((02)) 
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where C is the Euler constant, and # k is the angular de- 
pendence of the superconducting order parameter [4"5] , 
Eqn. 10 gives the transition temperature in paramag- 



netic, ferromagnetic and spiral phases provided \ and S 
are calculated in the appropriate backgrounds. We note 
that both the fluctuation-induced pairing and field renor- 
malization occur automatically in this approach. The 
field renormalization/self-energy term written here in- 
cludes the effect of the regularization described above 
Eqn. [3] Because of the intervention of phase reconstruc- 
tion before the quantum critical point is reached, we do 
not require that the magnetic fluctuations take their crit- 
ical form [37H2]. 

The dependence of E and x upon the background mag- 
netic order through the mean-field dispersion, e k (M, q) = 
e k + 0y/ (gM) 2 + (Ok ■ q) 2 , induces cross terms between 
the superconducting and magnetic order parameters and 
their gradients. Indeed such cross terms also appear in 
the mean-field expansion of the bare Cooper susceptibil- 
ity and ultimately in the pre-factor of chemical potential 
in Eq. ( 10 ) . Terms arising from the fluctuation corrections 



are precisely the mode-mode coupling terms included by 
Kirkpatrick et al. previously [4"0] . 

In practice, we find that although the gradient terms 
are crucial in driving the pair density wave state, ex- 
panding the susceptibility x an d self-energy £ in the spi- 
ral wave vector leads to sub-leading corrections that do 
not change the superconducting transition temperature 
appreciably. At small and large magnetization, we find: 
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PHASE DIAGRAM 

The phase diagram is shown in Fig. [l] With a first or- 
der transition and spin spiral emerging at low tempera- 
tures, the phase diagram is a refinement of that presented 
in Ref. [35] , the one notable difference being the order of 
the transition from the paramagnet to the spin spiral. In 
Ref. |35j only the lowest order term in the free energy ex- 
pansion was calculated, implicitly assuming a second or- 
der paramagnet-spiral transition. Here we have extended 
the free energy expansion to order M 6 , which shows that 
the paramagnet-spiral transition is in fact first order |20j . 
Though this phase diagram is consistent with variational 
Monte Carlo results [35], the spiral phase here occupies 
a smaller range of interaction strengths due to the ab- 
sence of sub-leading divergences. However the topology 
of the phase diagram is identical so we proceed with the 
amenable Ginzburg Landau formalism. As the supercon- 
ductivity is weak it may be laid over the magnetic phase 
diagram without appreciable feedback. 

Magnetic fluctuations drive Cooper pairing, so that 
even at M = a superconducting phase is formed. When 
the superconductivity overlaps regions of the phase di- 
agram supporting magnetic order, there is a large in- 
crease in the transition temperature because of the mode- 
mode coupling effects described above and in Ref. |40j . 
The transition temperature tails off exponentially to zero 
when the magnetization has depleted the number of pair- 
ing (spin-down) electrons to zero thus sending the sus- 
ceptibility term to zero. In Fig. [T] the superconduct- 
ing transition temperature shows a first order jump upon 
passing into the magnetically ordered phase, tracking the 
local magnetization. Where the superconductivity over- 
laps the spiral phase, its d-vector is locked to the local 
magnetic order forming a new phase, the pair density 
wave. 



DISCUSSION 

The minimal Hamiltonian Eq. - an electron gas with 
contact repulsion - delivers a remarkably rich phase di- 
agram when all possible fluctuations are accounted for. 
Fluctuations drive the itinerant ferromagnet near to sev- 
eral possible instabilities. Transverse fluctuations can 
drive a first order ferromagnetic transition or a spin spi- 
ral phase. Longitudinal magnetic fluctuations drive a 
p-wave superconducting instability. These instabilities 
overlap forming a pair density wave. Comparison with 
variational Monte Carlo calculations [35] shows that the 
quantum order-by-disorder approach, as well as reveal- 
ing the underlying physics in a transparent manner, can 
accurately determine the phase boundaries. 

However, the model may require modification in more 
realistic settings. Firstly, unlike the local interaction used 
here, the screened Coulomb interaction has finite range 
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and is sampled at different momenta by mean-field and 
fluctuation contributions to the free energy. The tricrit- 
ical point - the lowest temperature at which a continu- 
ous ferromagnetic to paramagnet transition is observed 
- is commonly seen at temperatures ~ 0.02Xf [T7] that 
are significantly lower than the prediction of the mini- 
mal Hamiltonian ~ 0.3Tp. However, it has been shown 
both with variational Monte Carlo and analytically that 
a finite ranged screened Coulomb interaction reduces 
the tricritical temperature found in a quantum order-by- 
disorder calculation to the realistic experimental range 
of ~ 0.02Tp 05]. Secondly, real band structures, with 
deviations from quadratic electron dispersions, shift the 
subtle balance between the near multi-critical phases [20] . 
Finally, spin-orbit induced magnetic anisotropy can sup- 
press transverse magnetic fluctuations relative to the lon- 
gitudinal fluctuations so changing the balance between 
magnetic and superconducting instabilities [38] , and may 
ultimately invalidate the neglect of back-reaction of su- 
perconducting order upon the magnetic order. 

The necessity of these modifications to the model 
Hamiltonian makes direct comparison of our results with 
experiment delicate. The types of order predicted in 
our model calculations have all been seen in (or sug- 
gested by) experiment. For example p-wave supercon- 
ducting is found in UGe 2 [3TH23], URhGe [23], and 
UCoGe [25]. In the magnetic state UGe2 displays su- 
perconductivity below T c sa 500mK, and no supercon- 
ductivity has been observed on the paramagnetic side. 
Theory predicts that the peak superconducting temper- 
ature has a ratio of 200 : 1 between the ferromagnetic 
and paramagnetic sides, meaning that on the param- 
agnetic side UGe2 would have a transition temperature 
at T c ss 3mK that is too small to measure, thus main- 
taining consistency with experiment [40 . Evidence for 
spin spiral behavior, though less direct, is found in Ce- 
FePO [Hi, NbFe 2 [J3], Yb(Rh .73Coo.27)2Si 2 [JJ] and 
possibly Sr 3 Ru207 [T71[TS]. Although evidence for the 
coexistence of superconductivity and spiral magnetism 
has not yet been found, it could be measured through a 
drop of electrical resistance coupled with the emergence 
of spatial texture in a neutron scattering study. Never- 
theless, we take the occurrence of these ingredients as 
reason for positive expectation. 

The quantum order-by-disorder framework provides a 
unifying and physical reformulation of diagrammatic ap- 
proaches to particle-hole and particle-particle instabili- 
ties. The Fermi surface provides boundary conditions for 
fluctuations through Pauli exclusion. Modification of the 
Fermi surface by the introduction of new order affects 
the spectrum of fluctuations, and can be self-consistently 
stabilized by them. We have shown that this drives the 
formation of a new phase, the pair density wave. This 
new formalism guides future calculations in a direction 
that may be fruitful in a variety of situations. 
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